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ABSTRACT 

This work aims to constrain the physical nature of umbral dots (UDs) using high-resolution spec- 
tropolarimetry. Full Stokes spectra recorded by the spectropolarimeter on Hinode of 51 UDs in a 
sunspot close to the disk center are analyzed. The height dependence of the temperature, magnetic 
field vector, and line-of-sight velocity across each UD is obtained from an inversion of the Stokes vec- 
tors of the two Fe I lines at 630 nm. No difference is found at higher altitudes (— 3 < log(r5oo) ^ — 2) 
between the UDs and the diffuse umbral background. Below that level the difference rapidly in- 
creases, so that at the continuum formation level (log(r5oo) = 0) we find on average a temperature 
enhancement of 570 K, a magnetic field weakening of 510 G, and upfiows of 800 m s~^ for peripheral 
UDs, whereas central UDs display an excess temperature of on average 550 K, a field weakening of 
480 G, and no significant upfiows. The results for, in particular, the peripheral UDs, including cuts of 
magnetic vector and velocity through them, look remarkably similar to the output of recent radiation 
MHD simulations. They strongly suggest that UDs are produced by convective upwellings. 
Subject headings: Sun: photosphere — Sun: sunspots — techniques: spectroscopic 



1. INTRODUCTION 

The energy transport immediately below the solar sur- 
face is mainly determined by convective processes that 
are visible as granulation patterns in white-light im- 
ages of the quiet photosphere. This convection is sup- 
pressed inside sunspot umbrae due to the strong ver- 
tical magn etic fi eld, but some form of magnetoconvec- 
tion ( Weissl l2^002f ) is needed to explain the observed um- 
bral brightnesses. Umbral fine structure such as light 
bridges or umbral dots, dot like bright features inside um- 
brae, may well be manifestations of magnetoconvection. 
Different models have been proposed to explain UDs, 
e.g., columns of field- free hot g:as in between a bundle 
of thin magnetic fiux ropes (jParked fTOTOl : IChoudhurvl 
[1986), or spatiall y modulated oscillations in a strong 
magnetic field (Weiss et al."l990) . Recent numerical sim- 
ulatio ns of three-dimensional radiative magnetoconvec- 
tion (ISchiissler fc VQg"Ierll2006l ) reveal convective plumes 
that penetrate through the solar surface and look very 
much like UDs. Although rece nt broadband images 
may h ave spat ially resolved UDs dSobotka fc Hansl meierl 
120051 : It. L. Riethmiiller et al.l [2008. in preparation), 
spectropolarimetry is needed to learn more about their 
physical nature. Previous spectroscopic observations led 
to heterogeneous results. Kneer (1973) found that UDs 
exhibit upfiows of 3 km s~^ and a 50% weaker magnetic 
field c ompa red to the nearby umb r a, whe reas Lites et al. 
(jl99lf ) and iTritschler fc Schmidtl (|1997D reporte d little 
field weakening. Finally, iSocas- Navarro et all (|2004f ) 
observed a weakening of 500 G and upfiows of a few 



obtained by the spectropolarimeter (SP) on the Hinode 
spacecraft. The usefulness of Hinode data for the study 
of UDs was demonstrated by Bharti et al.l (|2007) . who 
found that large UDs show dark lanes whose ex istence 
had been predicted bv ISchiissler fc Voglerl (|2006l ). 

2. OBSERVATIONS AND DATA REDUCTION 

The data em ployed here were acquired by the spec- 
tropolarimeter (|Lites et al] 1200 ll) of the Solar Optical 
?. ISuematsu et al.| [2 



100 m s" 



More details can be found in the reviews 



of urn bral fine structure bv ISolan"kil ([200 3) and Sobotka 
(j2006f ). One reason for the difference in results has been 
the infiuence of scattered light and variable seeing, which 
affect the different analyzed data sets to varying degrees. 
It therefore seems worthwhile to invert Stokes profiles 



Electronic address: [riethmueller;solanki;lagg]@mps. mpg.de 



Telescope (SOT. ISuematsu et al.ll2008[ ) onboard Hinode. 
They are composed of full Stokes spectra in the Fe I 
line pair around 6302 A and the nearby continuum of 
a sunspot of NOAA AR 10933 recorded from 12:43 to 
12:59 UT on 2007 January 5 using the 0.16''xl64'' slit. 
At this time the sunspot was located at a heliocentric an- 
gle of 4°, i.e. very close to disk center. The observations 
covered the spectral range from 6300.89 to 6303.26 A, 
with a sampling of 21 mA pixel" ^. The SP was operated 
in its normal map mode, i.e. both the sampling along 
the slit and the slit-scan sampling were 0.16'', so that 
the spatial resolution should be close to the diffraction 
limit of 1.22 X/D = 0.32''. The integration time per 
slit position was 4.8 s which reduced the noise level to 

10-3 

The data were corrected for dark current, fiat field, and 
instrumental polarization with the help of the SolarSoft 
package.^ A continuum intensity image (put together 
from the slit scan) of the chosen umbra is shown in Fig- 
ure [H Due to the large slit length we are always able to 
find a sufficiently extensive region of quiet Sun that is 
used to normalize intensities. 

3. DATA ANALYSIS 

To obtain atmospheric stratifications of temperature 
(T), magnetic field strength (5), and hne-of-sight ve- 
locity (vl os) we use th e inve rsion code SPINOR de- 
scribed bv lFrutiger et al.l (|2000l ). This code incorporates 

^ See http : //www . Imsal . com/solarsof t [ 
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Fig. 1. — Continuum intensity map of the sunspot NOAA 10933 
as observed by the Hinode SOT/SP on 2007 January 5. Hehocen- 
tric angle is ^ = 4°. Intensities are normahzed to the intensity level 
of the quiet photosphere 1^^. The white line at (4,4) Mm marks 
the cut through an umbral dot (UD) that is discussed in greater 
detail. 



the STOPRO routines (jSolankil [ToSTI ). which compute 
synthetic Stokes profiles of one or more lines upon in- 
put of their atomic data and one or more model atmo- 
spheres. Local thermodynamic equilibrium conditions 
are assumed and the Unno-Rachkovsky radiative transfer 
equations are solved. The inversions use an optical depth 
scale as the appropriate coordinate for radiative transfer 
problems. For reasons of comparability we use the op- 
tical depth at 500 nm (rsoo)- Starting with an initial 
guess model, the synthetic profiles were iteratively fitted 
to observed data us in^ response functions (RFs) and the 
merit function (|Ruiz Cobo fc del Toro Iniestal 1 19921 : 
lFrutigerll2QQQl ) is minimized. With the help of the RFs 
we find that the Fe I line pair at 6302 A is mainly formed 
within the log(r5oo) interval [—3, 0], which corresponds to 
a height range of about 400 km under hydrostatic equi- 
librium conditions in the umbra. The free parameters 
are defined at the four nodes —3, —2, —1, and of the 
log(r5oo) grid. The atmospheric stratification is then in- 
terpolated using splines onto a 10 times finer log(r5oo) 
grid. 

The first step of our analysis is the wavelength calibra- 
tion required to determine line-of-sight (LOS) velocities. 
For every slit position we average the Stokes / profiles 
of all locations along the slit whose total polarization 
P = /(Q^ + + V'^Y^'^dX is negligible, since those lo- 
cations are assumed to represent the quiet Sun. This 
mean / profile is used to fit Voigt profiles to the two 
Fe I lines from which the line center wavelengths are de- 
ter mined. The convective blueshift of 140 m s~^ (see 
iMartmez Fillet et al.l Il997l : iDravins et all Il98lf ) is then 
removed. 

The next step is to find an appropriate model atmo- 
sphere. Since we are interested in the atmospheric strati- 
fication of temperature, magnetic field strength, and LOS 
velocity within a UD, these three atmospheric parame- 
ters are assumed to be height dependent, whereas field 
inclination and azimuth angle, micro- turbulence, and 
macro-turbulence are assumed to be height independent. 
We experimented intensively with adding a second model 
component to represent the stray light, but the inversion 
results did not improve significantly, confirming the al- 




FlG. 2. — Stokes /, V ^ Q and U profiles from the center of the 
UD marked in Figure [T] Red lines are the measured, blue lines the 
best-fit profiles, i.e. the inversion result. The bottom parts of each 
panel show the difference between the two on an expanded scale. 

most negligible stray light in the SP. Therefore, in the 
interests of a robust inversion, we forbore from adding a 
stray light component, thus reducing the number of free 
parameters. 

Lastly, we have to find initial guesses for all free pa- 
rameters. We use an initial temperature stratifica tion 
according to the umbral core model L of Maltbv et al.l 
(Il986l ) and assume a vertical magnetic field of 2000 G and 
zero LOS velocity at all heights. Initial guesses for micro- 
turbulence and macro-turbulence are 0.1 and 2 km s~^, 
respectively. Other initial guesses gave very similar re- 
sults, except for a limited number of outliers. For these, 
repeating the inversion with an initial guess close to the 
final result of one of the neighboring pixels returned val- 
ues consistent with those obtained for the other pixels. 

4. INVERSION RESULTS 

We analyzed a total of 51 UDs, which were identi- 
fied by applying the mu l tilevel tracking (MLT) algorith m 
(jBovelet fc Wiehrll2001l : lT. L. Riethmiiller et al.ll2008l . in 
preparation). For each UD the location of its core was 
identified, a cut was made through it, reaching to the 
neighboring diffuse background (DB), and the profiles 
from all the pixels along this cut were inverted. We first 
discuss the results for the UD marked in Figure [H cho- 
sen because of its brightness, which leads to particularly 
small error bars. A comparison of the measured profiles 
with the best-fit profiles resulted from the inversion can 
be seen in Figure [2] for the UD and in Figure [3] for the 
DB selected as the location of lowest continuum inten- 
sity in a 1.4 X 1.4 Mm^ environment of the UD center. 
Due to the low signal in the dark background the mea- 
sured DB profiles are much noisier than the UD center's 
profiles, but in general, the Stokes spectra can be fitted 
remarkably well. 

The stratification of the retrieved atmospheric param- 
eters T, vloSi cind B in the center of the UD and in 
the DB are plotted in Figure [H In the upper photo- 
sphere (—3 < log(r5oo) < —2) the error bars overlap; 
i.e. we find little significant difference between UD and 
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Fig. 3. — The same as Figure [2] but for Stokes /, V, Q and U 
profiles of the diffuse background near the UD. 
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Fig. 4. — Atmospheric stratification obtained from the Stokes 
profiles at the location of the UD's center (red lines) and from 
the Stokes profiles of the diffuse background near the UD (blue 
lines). The formal errors of the inversion at the used optical depth 
nodes are indicated by bars. Negative LOS velocity values indicate 
upfiows. 

DB. In the deeper photosphere, however, the inversions 
return strongly different stratifications. Thus, the UD 
temperature is higher than the DB temperature, consis- 
tent with the intensity enhancement of the UD in the 
continuum map. The LOS velocity (which is identical 
to the vertical velocity due to the small heliocentric an- 
gle) exhibits strong upfiows in the UD center, whereas 
the DB is nearly at rest. The magnetic field strength is 
roughly 2 kG for the heights — 3 < log(r5oo) ^ — 1- Be- 
low log(r5oo) = — 1 the field strength of the UD decreases 
strongly with depth, whereas the field strength of the DB 
increases moderately. 

The vertical cuts of magnetic field strength and LOS 
velocity through 13 pixels lying along the white line in 
Figure [1] are shown in Figure [5l Jumps from one pixel to 
the next were smoothed through interpolation. There 
is clear evidence for a localized decrease in UD field 
strength in the low photosphere, co-located with an up- 
flow that extends higher into the atmosphere and a weak 
downflow on at least one side. The magnetic fields are 
4° more inclined in the UD than they are in the DB 
around t he UD. Fig ure [5] looks remarkably like Fig. 2 of 
ISchiissler fc Voglerl ([2006 ). in spite of the fact that Fig- 
ure [5] is plotted on an optical depth scale in the vertical 
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Fig. 5. — Vertical cut through the UD marked in Figure ^ in 
the direction indicated by the white line. Colors of the top panel 
indicate magnetic field strength. The bottom panel shows LOS 
velocity. Negative velocities are upfiows. 

direction and is thus distorted by an unknown amount 
relative to a corresponding figure on a geometrical scale. 

Next we discuss all 51 analyzed UDs. In the litera- 
ture we often find a separation into two UD regimes. 
For example, iGrossmann-Doerth et al.l (|l986l ) differen- 
tiate between peripheral UDs (PUDs) and central UDs 
(CUDs), i.e. between UDs that are born close to the 
umbra-penumbra boundary and UDs that are born deep 
in the umbra. We follow this distinction and plot the ob- 
tained stratifications of the 30 PUDs (distance to umbra- 
penumbra boundary less than 2000 km) in the top panels 
of Figureini while the remaining 21 CUDs are represented 
in the bottom panels of Figure [6l The results largely 
mirror those obtained for the UD discussed above. In 
the upper atmosphere UDs center and DB do not differ 
in their mean values of T, vloSi ^ind B. On average, 
the CUDs are about 150 K cooler than the PUDs in the 
upper atmosphere, just as the DB around the CUDs is 
cooler than the DB around the PUDs. At log(r5oo) = 
we find that PUDs are 570 K hotter than the local DB 
and CUDs are 550 K hotter than the DB in their vicin- 
ity. The magnetic field strength at log(r5oo) = is weak- 
ened by about 510 G for PUDs and 480 G for CUDs, 
whereas only PUDs exhibit significant upfiows of about 
800 m s~^. The mean LOS velocity shows no difference 
between CUD centers and DB. In order to make sure that 
an upflow is not being missed due to the lower S/N ratio 
of the CUD Stokes profiles, we have also averaged the 
Stokes profiles of all the CUDs. An inversion of there av- 
eraged Stokes profiles gave a result that agrees with the 
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Fig. 6. — Atmospheric stratifications of peripheral umbral dots 
(top 3 panels) and central umbral dots (bottom 3 panels). The red 
lines show the stratification at the location of the UD's center and 
the blue lines correspond to the nearby diffuse background. The 
green line is the weighted average of all red lines and the yellow 
line is the weighted average of all blue lines, where we used the 
reciprocal error bars as weighting factors. 

averaged stratifications {green line) in the bottom pan- 
els of Figure [6] within the error bars. This suggests that 
any upflow velocity in CUDs is mostly restricted to layers 
below the surface or is too concentrated or too weak to 
be detected by the inversions. Finally, we find that the 
magnetic field of the PUDs is on average 4° more hori- 
zontal than for their DB. We see no inclination difference 
for CUDs. 

5. DISCUSSION 

We identified 30 peripheral and 21 central umbral dots 
in Hinode spectropolarimetric data of a sunspot within 
4° of disk center. With the help of Stokes profile in- 
versions of the Fe I lines at 630 nm we determined the 
stratifications of temperature, magnetic field strength. 



an d LOS velocity. Th e pres ent work differs from that 
of iSocas-Navarro et al.l (|2004l ) in the superior quality of 
the employed data with twice the spatial resolution and 
practically no scattered light. This allows a detailed de- 
termination of the atmospheric stratification. The higher 
spatial resolution of the Hinode SP data also allows us 
to, for the first time, reconstruct both the horizontal and 
the vertical structure of UDs. We also extended the anal- 
ysis to a more numerous statistical ensemble of 51 UDs. 

Vertical cuts through UDs provide a remarkable 
confirmation of th e resu lts of MHD simulations of 
ISchiissler Voglerl (|2006f ): both show that UDs dif- 
fer from their surroundings mainly in the lowest visi- 
ble layers, where the temperature is enhanced and the 
magnetic field is weakened. We found a temperature 
enhancement of 550 K and a magnetic field reduction 
of about 500 G (at optical depth unity). In addition, 
PUDs display upfiow velocities of 800 m s~^ on average, 
again in good agreement with the simulations. There 
are also some differences between our results and those of 
ISchiissler & Vogler (2006). Thus, according to our inver- 
sions the magnetic field strength of the DB is somewhat 
depth dependent. This was not the case for the MHD 
simulations due to the used periodic boundary condi- 
tions. Furthermore, although some of the UDs display a 
weak downfiow bounding the strong central upfiow (see 
Figure [5]) , these are neither as narrow nor as strong as 
the downfiows at the end s of dark lanes as reported by 
ISchiissler fc Voglerl (|2006l ). probably due to the limited 
spatial resolution of our data. We may also be missing 
some of the narrow downfiows by considering only single 
c uts across individual UDs. 

I So cas- Navarro et al] (|2004l ) reported 10° more inclined 
magnetic fields in PUDs. This result is qualitatively 
confirmed by our work; we find an inclination increase 
of 4° for PUDs but no increase for CUDs, which can 
be assumed as a further hint that the main part of the 
CUD structure is below the surface. These results can 
be interpreted in terms of the strong DB fields expand- 
ing with height and closing over the UD, as proposed by 
■Socas-Navarro et al., (2004 ). 



REFERENCES 



Bharti, L., Joshi, C. & Jaaffrey, S. N. A. 2007, ApJ, 669, L57 
Bovelet, B., & Wiehr, E. 2001, Solar Phys., 201, 13 
Choudhury, A. R. 1986, ApJ, 302, 809 

Dravins, D., Lindegren, L., & Nordlund, A. 1981, A&A, 96, 345 
Frutiger, C. 2000, Inversion of Zeeman Spht Stokes Profiles: 
Application to solar and stellar surface structures, Ph.D. Thesis, 
Institute of Astronomy, ETH Zurich, No. 13896 
Frutiger, C, Solanki, S. K., Fligge, M., & Bruls, J. H. M. J. 2000, 
A&A, 358, 1109 

Grossmann-Doerth, U., Schmidt, W., & Schroter, E. H. 1986, A&A, 
156, 347 

Kneer, F. 1973, Solar Phys., 28, 361 

Lites, B. W., Bida, T. A., Johannesson, A., Sz Scharmer, G. B. 

1991, ApJ, 373, 683 
Lites, B. W., Elmore, D. F., & Streander, K. V. 2001, in ASP 

Conf. Ser. 236, Advanced Solar Polarimetry, ed. M. Sigwarth 

(San Francisco: ASP), 33 
Maltby, P., Avrett, E. H., Carlsson, M., Kjeldseth-Moe, O., Kurucz, 

R. L., et al. 1986, ApJ, 306, 284 
Martmez Fillet, V., Lites, B. W. &; Skumanich, A. 1997, ApJ, 474, 

810 

Parker, E. N. 1979, ApJ, 234, 333 



Riethmiiller, T. L., Solanki, S. K., Zakharov, V. &; Gandorfer, A. 

2008, A&A, 492, 233 
Ruiz Cobo, B., &; del Toro Iniesta, J. C. 1992, ApJ, 398, 375 
Schiissler, M., &; Vogler, M. 2006, ApJ, 641, L73 
Sobotka, M. 2006, Dissertation for Doctor Scientiarum, Acad. Sci. 

Czech Republic 
Sobotka, M., &; Hanslmeier, A. 2005, A&A, 442, 323 
Socas-Navarro, H., Martmez Fillet, V., Sobotka, M., & Vazquez, 

M. 2004 ApJ, 614, 448 
Solanki, S. K. 1987, The Photospheric Layers of Solar Magnetic 

Fluxtubes, Ph.D. Thesis, Institute of Astronomy, ETH Ziirich, 

No. 8309 

Solanki, S. K. 2003, A&A Rev., 11, 153 

Suematsu, Y., Tsuneta, S., Ichimoto, K., Shimizu, T., Otsubo, M., 

et al. 2008, Solar Phys., in press 
Tritschler, A., &; Schmidt, W. 1997, A&A, 321, 643 
Weiss, N. O. 2002, Astron. Nachr., 323, 371 

Weiss, N. O., Brownjohn, D. P., Hurlburt, N. E., & Proctor, M. R. 
E. 1990, Mon. Not. R. astr. Soc, 245, 434 



